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I. INTRODUCTION 



the matter and A contribution to the total energy den- 
sity. As a response to these problems, much interest has 
been devoted to models with dynamical vacuum energy, 
dubbed quintessence [lfj . These models typically involve 
a scalar field rolling down its self interaction potential 
thus allowing the vacuum energy to become dominant 
only recently (see ^3, El for good reviews) . Although 
quintessence by a scalar field is the most studied candi- 
date for dark energy, it generally does not avoid ad hoc 
fine tuning to solve the coincidence problem. Moreover, 
it is not clear where this scalar field comes from and how 
to choose the self interaction potential. 

On the other hand, it is worth noting that, despite 
the broad interest in dark matter and dark energy, their 
physical properties are still poorly understood at a fun- 
damental level and, indeed, it has never been shown that 
they are in fact two different ingredients. This obser- 
vation motivated the great interest recently devoted to 
a completely different approach to quintessence. Rather 
than fine tuning a scalar field potential, it is also possible 
to explain the acceleration of the universe by introducing 
a cosmic fluid with an exotic equation of state causing it 
to act like dark matter at high density and dark energy at 
low density. An attractive feature of these models is that 
they can explain both dark energy and dark matter with 
a single component (thus automatically solving the coin- 
cidence problem) and have therefore been referred to as 
unified dark energy (UDE) or unified dark matter (UDM). 
Some interesting examples are the generalized Chaplygin 
gas Il3j . the tachyonic field Q, the condensate cosmol- 
ogy jl5j| and the Hobbit model 0] . It is worth noting, 
however, that such models are seriously affected by prob- 
lems with structure formation 01 so that some work is 

•Electronic address: capozziello@sa.infn.it, winny@na.infn.it, francavigl^gPweisieditiefore they Can be considered as reliable al- 
1 It is common in literature to make no distinction between the ternatives to dark energy 
concordance and the ACDM model even if, strictly speaking, in 

Actually, there is still a different way to face the prob- 



The Hubble diagram of type la supernovae (here- 
after SNela) P, , the anisotropy spectrum of the cos- 
mic microwave background radiation (hereafter CMBR) 
[J 0] the matter power spectrum determined by the 
large scale distribution of galaxies [|J and by the data 
on the Lya clouds [8j are all convincing evidences in 
favour of a new picture of the universe, unexpected only 
few years ago. According to this nowadays standard sce- 
nario, the universe is flat and undergoing an accelerated 
expansion driven by a mysterious fluid with negative 
pressure nearly homogeneously distributed and making 
up to ~ 70% of the energy content. This exotic compo- 
nent is what is called dark energy, while the model we 
have just depicted is usually referred to as the concor- 
dance model. 

Even if strongly supported by the bulk of the available 
astrophysical data, this new picture is not free of prob- 
lems. Actually, while it is clear how dark energy works, 
its nature remains an unsolved problem. The simplest 
explanation claims for the cosmological constant A thus 
leading to the so called ACDM model 1 9]. Although 
being the best fit to most of the available astrophysical 
data 0, , the ACDM model is also plagued by many 
problems on different scales. If interpreted as vacuum 
energy, A is up to 120 orders of magnitudes smaller than 
the predicted value. Furthermore, one should also solve 
the coincidence problem, i.e. the nearly equivalence of 



the concordance model the dark energy may also be provided by 
a different mechanism. 



lem of cosmic acceleration. As stressed in Lue et al. Ilf 
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it is possible that the observed acceleration is not the 
manifestation of another ingredient in the cosmic pie, 
but rather the first signal of a breakdown of our under- 
standing of the laws of gravitation. From this point of 
view, it is thus tempting to modify the Friedmann equa- 
tions to see whether it is possible to fit the astrophysical 
data with a model comprising only the standard mat- 
ter. Interesting examples of this kind are the Cardassian 
expansion and the DGP gravity pol |. 

In this same framework, there is also the attractive 
possibility to consider the Einsteinian general relativity 
as a particular case of a more fundamental theory. This 
is the underlyi ng philosophy of what are referred to as 
f(R) theories [Mill El 111 El E| In this case, the 
Friedmann equations have to be given away in favour of a 
modified set of cosmological equations that are obtained 
by varying a generalized gravity Lagrangian where the 
scalar curvature R has been replaced by a generic func- 
tion f(R). The usual general relativity is recovered in 
the limit f(R) = R, while completely different results 
may be obtained for other choices of f(R)- While in the 
weak field limit the theory should give the usual new- 
tonian gravity, at cosmological scales there is an almost 
complete freedom in the choice of f(R) thus leaving open 
the way to a wide range of models. 

The key point of f(R) theories is the presence of modi- 
fied Friedmann equations obtained by varying the gener- 
alized Lagrangian. However, here lies also the main prob- 
lem of this approach since it is not clear how the varia- 
tion has to be performed. Actually, once the Robertson - 
Walker metric has been assumed, the equations govern- 
ing the dynamics of the universe are different depending 
on whether one varies with respect to the metric only or 
with respect to the metric components and the connec- 
tions. It is usual to refer to these two possibilities as the 
metric and the Palatini approach respectively. The two 
methods give the same result only in the case f(R) = R, 
while lead to significantly different dynamical equations 
for every other choice of f(R) (see |!l El El (HI and 
references therein). 

It is worth noting f(R) theories were initially investi- 
gated using the metric approach [^lElES] • Even if some 
interesting and successful results have been obtained [3(J , 
this way to f(R) theories is plagued by serious mathe- 
matical difficulties. Actually, even for the simplest f(R), 
the metric approach leads to a fourth order nonlinear dif- 
ferential equation for the scale factor that is impossible 
to solve analytically and is affected by several problems 
that greatly complicate the search for numerical solu- 
tions. Moreover, some doubts have been cast on the con- 
sistency among the weak field limit of the theory and the 
newtonian gravity as tested at the Solar system scale jU 
even if some interesting different results have also been 
obtained |3^|. 

On the other hand, theoretical considerations about 
the stability of the equations and the newtonian limit 
argue in favor of the Palatini approach to f(R) theo- 
ries. Moreover, the dynamics of the universe may be 



analytically determined from the cosmological equations 
obtained with this method for some interesting cases. To 
this aim, a clear mathematical machinery has been pre- 
sented in Ref. EU (hereafter ABF04) that allows to de- 
termine analytic expressions for the Hubble parameter as 
function of the redshift. As we will see later, this is all 
what is needed to test a given cosmological model. 

The Palatini approach to f(R) theories has been 
widely studied in literature El il El El El El and 
the dynamics of the cosmological models obtained by ap- 
plying this method to different choices of f(R) has been 
investigated in detail. Here we adopt an observational 
point of view on the Palatini approach. Assuming that 
this is the correct way to treat f(R) theories, we in- 
vestigate the viability of two classes of models obtained 
by two popular choices for f(R), namely the power law 
f(R) = (3R n and the logarithmic f{R) =abxR. To this 
aim, we compare the model predictions against the SNela 
Hubble diagram and the data on the gas mass fraction 
in relaxed galaxy clusters. This analysis will allow us to 
constrain the model parameters and to see whether f(R) 
theories are indeed reliable alternatives to dark energy. 
Moreover, this will be an observational validation of the 
theoretically motivated Palatini approach. 

The paper is organized as follows. Sect. II details the 
method we employ to constrain the models and present 
the dataset we will use. The two classes of models we 
consider are briefly discussed in Sect. Ill where we also 
individuate the parameters that are better suited to both 
assign the model and be constrained by the data. A 
detailed discussion of the results is the subject of Sect. IV, 
while we summarize and conclude in Sect. V. 



II. CONSTRAINING A MODEL 

Considered for a long time a purely theoretical science, 
cosmology has today entered the realm of observations 
since it is now possible to test cosmological predictions 
against a meaningful set of astrophysical data. Much 
attention, in this sense, has been devoted to standard 
candles, i.e. astrophysical objects whose absolute magni- 
tude M is known (or may be exactly predicted) a priori 
so that a measurement of its apparent magnitude m im- 
mediately gives the distance modulus /x = m — M . The 
distance to the object is then estimated as : 



H{z) = 51og£ L (z) + 25 (1) 

with Dl(z) the luminosity distance (in Mpc) and z the 
redshift of the object. The relation between /x and z is 
what is referred to as Hubble diagram and is an open win- 
dow on the cosmography of the universe. Furthermore, 
the Hubble diagram is a powerful cosmological test since 
the luminosity distance is determined by the expansion 
rate as : 
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D L (z) = — (1 
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with -E(-z) = H(z)/Hq, H — d/a the Hubble parameter 
and a(t) the scale factor. Note that an overdot means 
differentiation with respect to cosmic time, while an un- 
derscript denotes the present day value of a quantity. 

Being the Hubble diagram related to the luminosity 
distance and being Dl determined by the expansion rate 
H[z), it is clear why it may be used as an efficient tool 
to test cosmological models and constrain their param- 
eters. To this aim, however, it is mandatory that the 
relation /i = ix{z) is measured up to high enough rcdshift 
since, for low z, Dl reduces to a linear function of the 
redshift (thus recovering the Hubble law) whatever the 
background cosmological model is. This necessity claims 
for standard candles that are bright enough to be visible 
at such high redshift that the Hubble diagram may dis- 
criminate among different rival theories. SNela are, up 
to now, the objects that best match these requirements. 
It is thus not surprising that the first evidences of an 
accelerating universe came from the SNela Hubble dia- 
gram Q] and dedicated survey (like the SNAP satellite 
[23) have been planned in order to increase the number 
of SNela observed and the redshift range probed. 

The most updated and reliable compilation of SNela is 
the Gold dataset recently relased by Riess et al. • The 
authors have compiled a catalog containing 157 SNela 
with z in the range (0.01,1.70) and visual absorption 
Ay < 0.5. The distance modulus of each object has been 
evaluated by using a set of calibrated methods so that the 
sample is homogenous in the sense that all the SNela have 
been re-analyzed using the same technique in such a way 
that the resulting Hubble diagram is indeed reliable and 
accurate. Given a cosmological model assigned by a set 
of parameters p = (jp\, . . . ,p n ), the luminosity distance 
may be evaluated with Eq.© and the predicted Hubble 
diagram contrasted with the observed SNela one. Con- 
straints on the model parameters may then be extracted 
by mean of a \ 2 ~ based analysis defining the \ 2 as : 
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where <7i is the error on the distance modulus at rcdshift 
Zi and the sum is over the NsNeia SNela observed. It 
is worth stressing that the uncertainty on each measure- 
ment also takes into account the error on the redshifit 
and are not gaussian distributed. As a consequence, the 
reduced \ 2 (i- e -i XsNeia divided by the number of degrees 
of freedom) for the best fit model is not forced to be close 
to unity. Nonetheless, different models may still be com- 
pared on the basis of the x 2 value : the lower is XsNeia> 
the better the model fits the SNela Hubble diagram. 

The method outlined above is a simple and quite ef- 
ficient way to test whether a given model is a viable 



candidate to describe the late time evolution of the uni- 
verse. Nonetheless, it is affected by some degeneracies 
that could be only partially broken by increasing the 
sample size and extending the redshift range probed. 
A straightforward example may help in elucidating this 
point. Let us consider the flat concordance cosmological 
model with matter and cosmological constant. It is : 

e 2 (z) = n M (i + z ) 3 + (i-n M ) 

so that x'sNeia w * u on ly depend on the Hubble constant 
Hq and the matter density parameter Qm- Actually, we 
could split the matter term in a baryonic and a non bary- 
onic part denoting with fib the baryon density parameter. 
Since both baryons and non baryonic dark matter scales 
as (1 + z) 3 , E(z) and thus the luminosity distance will de- 
pend only on the total matter density parameter and we 
could never constrain flf, by fitting the SNela Hubble dia- 
gram. Similar degeneracies could also happen with other 
cosmological models thus stressing the need for comple- 
mentary probes to be combined with the SNela data. 

To this aim, we consider a recently proposed test based 
on the gas mass fraction in galaxy clusters. We briefly 
outline here the method referring the interested reader 
to the literature for further details [3 HH • Both theo- 
retical arguments and numerical simulations predict that 
the baryonic mass fraction in the largest relaxed galaxy 
clusters should be invariant with the redshift (see, e.g., 
Ref. [3?! ). However, this will only appear to be the case 
when the reference cosmology in making the baryonic 
mass fraction measurements matches the true underly- 
ing cosmology. From the observational point of view, it is 
worth noting that the baryonic content in galaxy clusters 
is dominated by the hot X - ray emitting intra-cluster gas 
so that what is actually measured is the gas mass fraction 
f gas and it is this quantity that should be invariant with 
the redshift within the caveat quoted above. Moreover, it 
is expected that the baryonic fraction in clusters equals 
the universal ratio Qf,/Q,M so that f gas should indeed be 
given by bxQb/^M where the multiplicative factor b is 
motivated by simulations that suggest that the gas frac- 
tion is slightly lower than the universal ratio because of 
processes that convert part of the gas into stars or eject 
it outside the cluster itself. 

Following Ref. [3(| (hereafter A04), we adopt the 
SCDM model (i.e., a flat universe with Qm = 1 and 
h — 0.5, being h the Hubble constant in units of 
100 km s _1 Mpc -1 ) as reference cosmology in making 
the measurements so that the theoretical expectation for 
the apparent variation of f gas with the redshift is j3|| : 
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where D^ CDM and D™ od is the angular diameter dis- 
tance for the SCDM and the model to be tested respec- 
tively. Da(z) may be evaluated from the luminosity dis- 
tance Dl{z) as : 
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D A (z) = (l + z)- 2 D L (z) 



(5) 



with Dl(z) given by Eq.(J2J) above. 

A04 have extensively analyzed the set of simulations in 
Ref. jjtfl to get b = 0.824±0.089. In our analysis below, 
we will set b = 0.824 in order to not increase the num- 
ber of parameters to be constrained. Actually, we have 
checked that, for values in the la range quoted above, 
the main results are independent on b. It is worth not- 
ing that, while the angular diameter distance depends 
on E(z) and thus on h and VLm, the prefactor in Eq.QJ 
makes f gas explicitly depending on S} /£Im so that a di- 
rect estimate of fl is (in principle) possible. Actually, 
we will see later that, for the models we will consider, 
the quantity that is constrained by the data is the ratio 
51b/£ljv/ rather than Q itself. 

To simultaneously take into account both the fit to the 
SNela Hubble diagram and the test on the f gas data, it is 
convenient to perform a likelihood analysis defining the 
following likelihood function : 



£(p) oc exp 



x 2 (p) 



(6) 



with : 



SNela Hubble diagram, the CMBR anisotropy spectrum 
measured by WMAP and the matter power spectrum ex- 
tracted from over 200000 galaxies observed by the SDSS 
collaboration. It is worth noting that, while our prior on 
h is the same as that used by many authors when apply- 



ing the f gas test 35, 36], it is common to put a second 
prior on Q rather than £1 /VIm- Actually, this choice 
is motivated by the peculiar features of the models we 
will consider that lead us to choose this unusual prior for 
reasons that will be clear later. 

With the definition © of the likelihood function, the 
best fit model parameters are those that maximize C(p). 
However, to constrain a given parameter pi, one resorts 
to the marginalized likelihood function defined as : 



£ Pi (Pi) oc dpi... / dpi^x / dp i+ i 



dp n £(p) (9) 



that is normalized at unity at maximum. The la con- 
fidence regions are determined by S\ 2 = X 2 ~ Xo = 1> 
while the condition 5x 2 = 4 delimited the 2a confidence 
regions. Here, Xo ^ s the value of the \ 2 f° r the best fit 
model. Projections of the likelihood function allow to 
show eventual correlations among the model parameters. 
In these two dimensional plots, the la and 2a regions are 
formally defined by A\ 2 = 2.30 and 6.17 respectively so 
that these contours are not necessarily equivalent to the 
same confidence level for single parameter estimates. 
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where we have defined : 
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being /^(zj) the measured gas fraction in a galaxy clus- 
ters at redshift Zi with an error a g i and the sum is over 
the Ngas clusters considered. In order to avoid possible 
systematic errors in the f gas measurement, it is desirable 
that the cluster is both highly luminous (so that the S/N 
ratio is high) and relaxed so that both merging processes 
and cooling flows are absent. A04 [3(| have recently re- 
leased a catalog comprising 26 large relaxed clusters with 
a precise measurement of both the gas mass fraction f gas 
and the redshift z (not presented in the quoted paper). 
We use these data to perform our likelihood analysis in 
the following. 

Note that, in Eq.JZJ, we have explicitly introduced two 
gaussian priors to better constrain the model parameters. 
First, there is a prior on the Hubble constant h deter- 
mined by the results of the HST Key project j^j from 
an accurate calibration of a set of different local distance 
estimators. Further, we impose a constraint on the ratio 
fift/Ojvf by considering the estimates of Q, h 2 and ilmh 2 
obtained by Tegmark et al. Q from a combined fit to the 



III. THE f(R) MODELS 

The observed cosmic acceleration is currently ex- 
plained by invoking the presence of a new fluid with 
negative pressure which smoothly fills the universe driv- 
ing its expansion. However, the nature and the nurture 
of this fluid are yet unknown so that other radically 
different proposals, such as unified dark energy models 
El III EE IS or brane world inspired theories El US, 
are still viable and worth exploring. 

A quite interesting and fascinating scenario predicts 
that standard matter is the only ingredient of the cos- 
mic pie as it is indeed observed, but the Einsteinian gen- 
eral relativity breaks down at the present small curvature 
scale. As a result, one should generalize the action as: 



A = J [y/gf(R) + 2 K L mat ] d 4 



with k — 8nG and L mat the matter Lagrangian. Vary- 
ing with respect to the metric components and adopting 
then the Robertson - Walker metric, one obtains modi- 
fied Friedmann equations that, by rearranging suitably 
the different terms, may still be formally written in the 
same way as the usual ones provided that a new fictitious 
component is added. For instance, the Hubble parameter 
is now given as : 



H 2 
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(10) 
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with p m the standard matter energy density and p curv 
the energy density of a curvature fluid whose density and 
pressure are given in terms of f(R) and its derivatives 
(see 0, |23j for details). Although intriguing, this ap- 
proach leads to a mathematically untractable problem. 
Indeed, it turns out that the scale factor a(t) should be 
obtained by solving a nonlinear fourth order differential 
equation. Not surprisingly, it is not possible to analyti- 
cally solve this equation even for the simplest choices of 
f{R). Moreover, some conceptual difficulties make it a 
daunting task to look for numerical solutions. 

An attractive way to escape these problems is to re- 
sort to the so called Palatini approach in which the field 
equations are obtained by varying with respect to both 
the metric components and the connections considered 
as independent variables. A consistency condition is then 
imposed to complement the system thus giving a set of 
first order differential equations for the scale factor a(t) 
and the scalar curvature R. The modified Friedmann 
equations are finally obtained by imposing that the met- 
ric is the Robertson- Walker one (see, e.g., 28] for a clear 
illustration of the procedure). 

The Palatini approach is physically well motivated and 
has the attractive feature that the Hubble parameter 
H(z), that is all what is needed for constraining the 
model, may be expressed analytically for some choices 
of the function f(R). It is thus quite interesting to con- 
strain the cosmological models obtained by applying the 
Palatini approach with two different choices of the func- 
tion f(R). The main characteristics of these models are 
briefly presented below. We follow Ref. (hereafter 
ABF04) which the interested reader is referred to for fur- 
ther details. 
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with r\ = p m (z = 0) the present day value of the mat- 
ter density and e = ±1 depending on n in such a way 
that both a(t) and H(z) are correctly defined. For the 
applications, it is better to use the following relation : 

KT] = 3£ImHq 

with VLm the usual matter density parameter. It is worth 
stressing that, even if we assume a flat model, £Im is not 
forced to be unity since the critical density for closure 
is now different from the usual value p c — 3Hq/8itG. 
The present day age of the universe may be obtained by 
evaluating Eq. l|12[l at the present day and then solving 
with respect to to thus obtaining : 



(14) 



Being the scale factor a power law function of the time, 
the deceleration parameter is constant and given as : 
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A. The Power law Lagrangian 

We first consider the class of Lagrangians that are lin- 
ear in an aribtrary power of the scalar curvature R : 



so that we may exclude all the Lagrangians with n < 3/2 
since they give rise to non accelerating models (go > 0). 

A nice feature of this model is that the dimensionless 
Hubble parameter is simply : 



W) = PR 1 ' 



(ii) 



with (3 7^ and n ^ 0, 2 real parameters to be con- 
strained. Note that (3 has the same units of R n so 
that f(R) is adimensional. This model has been already 
discussed by many authors 0, H^] using the stan- 
dard way of varying the Lagrangian. In particular, in 
Ref. [33 , some of us have also successfully tested a simpli- 
fied version of this model (with no matter term) against 
the SNela Hubble diagram. Moreover, this kind of La- 
grangian has also been investigated in the framework of 
the Palatini approach [24L |28| . It is thus particularly in- 
teresting to see whether the Palatini approach leads to 
results that are in agreement with the observed data. Us- 
ing the same notation as in ABF04, the scale factor a(t) 
and the Hubble parameter H(z) for a flat universe are 
given as : 



E 2 {z) = {l + zf /n (16) 
so that the luminosity distance turns out to be : 



Dl(z) 



2n 



Hn2n~3 



(1 + z) (1 + z) 



1 



(17) 



Both Dl and Da = (1 + z)~ 2 Dl depend only on the 
two parameters n and Hq so that fitting to the SNela 
Hubble diagram is unable to put any constraint neither 
on (3 or ft m ■ Adding the test on the f gas data described 
in the previous section partially alleviates this problem 
since f ga s{z) depends also on fib/ fin- It is then possi- 
ble to get an estimate of fijvf combining the constraint 
on Qb/rijM with an independent knowledge of from 
the measured abundance of light elements or primordial 
nucleosynthesis. Finally, the coupling parameter could 
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be derived inverting Eq. (|14f> with respect to (3 itself pro- 
vided that t has been somehow evaluated (possibly from 
a model independent method). 

As a general remark, let us observe that, without a 
knowledge of to, the parameter that can be constrained 
is Hm/P- Qualitatively, this could be explained by not- 
ing that all the tests we are considering are related to the 
cosmography of the universe. This is determined by the 
balance between the matter content and the exotic geo- 
metrical effects due to the replacement of R with f(R) in 
the gravity Lagrangian. Actually, this feature is common 
to all f(R) theories and could be expected since now ge- 
ometry plays the same role as the scalar field in the usual 
dark energy models. 



B. The Logarithmic Lagrangian 

Quantum effects in curved spacetimes may induce log- 
arithmic terms in the gravity Lagrangian 25] . It is thus 
interesting to consider the choice : 



f(R) = alni? 



(18) 



where the dimensions of a are such that f{R) is dimen- 
sionless 2 . This model is more complicated than the power 
law one so that, as a result, it is not possible to derive an 
analytical expression for the scale factor. However, the 
dimensionless Hubble parameter may still be expressed 
analytically as : 



E 2 {z) 
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(9/4)f» M #o 2 a~ 



(9/4)fi M -ff 2 a- 1 (l + ^) 3 



1 



x exp {(3/2)fWffoa _1 [(1 + z f - l] } -(19) 

The luminosity density is obtained inserting Eq. I jlfljl into 
the definition (J2J. There is not an analytic expression 
for Dl, but the integral is straightforward to evaluate 
numerically for a given value of flj^Hoa . As a con- 
sequence, the likelihood function for this model depend 
on the Hubble constant Hq, the ratio fib /fin between 
the baryonic and total matter density and the combined 
parameter QmHqCx - . It is worth stressing that, even 
if in principle possible, constraining separately the three 
parameters (flM,H ,a) is not correct since both Dl{z) 
and f g as(z) depend on a only through the combination 
f2jvfijga~ . Henceforth, it is this quantity that is con- 
strained by the data. Actually, this degeneracy may be 



broken by an independent estimate of fib that can be 
combined with the constraint on Q^/Qm to evaluate flj^ 
and then a from the constrained Q,mHqO. ~ l . Note that, 
without an estimate of Of, the only quantities estimated 
from the fit to the SNela Hubble diagram are Hq and 
flhiH^a^ 1 so that only the parameter fl m /a may be con- 
strained as a result of the above mentioned degeneracy 
between matter and geometry. 

There is no explicit analytic expression for the age of 
the universe so that one has to resort to numerical inte- 
gration of the following relation : 



in = 9.78 h~ 



d( 



(1 + C)#(C) 



(20) 



giving to expressed in Gyr. Le us remark that, while for 
power law Lagrangians to and fib are needed to break the 
degeneracy flu/ P, now fib and the likelihood analysis are 
sufficient to estimate both flM and a so that to may be 
used to check the results against an independent quantity. 

Another striking difference with the case of power law 
f(R) is the fact that the deceleration parameter is no 
longer constant. Even if we do not have an analytic ex- 
pression for a(t), we may still evaluate q as follows : 



-1 



1 + zdH 
H dz 



Inserting Ea. (|19|l into the above relation and evaluating 
the result at the present day (z = 0), we get : 



54 



108 



^) (21) 



l + 9fl M H £ a- 1 4 + 9fl M H { 

Ea. (l21|) shows that qo depends only on the parameter 
flMH^a' 1 that is therefore what determines whether the 
universe is today accelerating or decelerating. It is also 
worth noting that q{z) (not explicitly reported here for 
sake of shortness) changes sign during the evolution of 
the universe so that it is possible to estimate a tran- 
sition redshift zt as q(zr) = that only depends on 
flAiHo® 1 ■ ^ should be possible to estimate somewhat 
zt, this could give an independent check of the results. 
Actually, we will see that this is not possible since all 
the estimates of zt are model dependent. However, it is 
interesting to compare the transition redshift predicted 
for the logarithmic f(R) with that of other dark energy 
models. 



IV. RESULTS 



2 Note that, in literature, it is sometimes adopted the choice 
f(R) = alnPiR. We follow ABF04 and set ft = 1 with no 
loss of generality. 



We have applied the method described in Sect. II to 
investigate whether the cosmological models obtained by 
applying the Palatini approach to f(R) theories for the 
two choices in Eas. l|ll[l and i|18|) are in agreement with 
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both the SNela Hubble diagram and the data on the gas 
mass fraction in relaxed galaxy clusters. This also allows 
us to constrain the model parameters and compare the 
estimated values of some of them (as the Hubble constant 
h and the matter density Qm) with the recent results in 
literature in order to see whether they are reliable or not. 

A. f(R) = (5R n 

Let us first discuss the case of the power law La- 
grangian. The best fit parameters turn out to be : 

n = 2.25 , h = 0.641 , Q 6 /fi M = 0.181 (22) 

that gives the best fit curves shown in Figs.n an d[21 The 
agreement with the data (in particular, with the SNela 
Hubble diagram) is quite good which should be consid- 
ered a strong evidence in favor of the model. However, 
Fig. H shows that the model slightly overpredicts the dis- 
tance modulus for two highest redshift SNela, but, given 
the paucity of the data in this redshift range, the discrep- 
ancy is hardly significative. Should this trend be con- 
firmed by future data (observable, e.g., with the SNAP 
satellite mission that will detect SNela up to z ~ 2), we 
should exclude the choice (fTT|> for f(R). Actually, such 
a result could be expected since the deceleration param- 
eter is constant, while Riess et al. 2] claimed to have 
detected a change in the sign of q at a transition redshift 
zt ~ 0.5. We will return later to the problems connected 
with the result of Riess et al. that lead us to consider (at 
least) premature to deem as unreliable a model with a 
constant q. Therefore, we still retain f(R) theories with 
power law Lagrangian. 

It is interesting to look at the confidence contours 
in the projected two parameters space. Figs.|31 and 0] 
show the confidence regions for the parameters (n, h) and 
(n, Qb/£lM) respectively. It turns out that n is positively 
correlated with both h and flb/flm so that the higher 
is n, the higher is the expansion rate and the lower is 



the matter content fl^j- As a consequence, to fit the 
available data, models with steeper (higher n) power law 
Lagrangians should contain less matter which is a result 
disfavoring values of n much larger than our best fit. 

Using the method detailed in Sect. II, we have obtained 
the following constraint on the model parameters : 

f (2.06,2.46) at la 
nel ; (23) 

[ (1.91,2.61) at 2a 

( (0.637,0.648) at la 
he I ; (24) 

[ (0.633,0.654) at 2a 

( (0.177,0.185) at la 

n h /n M e { . (25) 

[ (0.173,0.189) at 2a 

The cosmological model originating from power law f(R) 
has been already considered by different authors in liter- 
ature under the metric approach to the variation of the 
Lagrangian pH I22I |23| . However, the lack of analytic so- 
lutions for the scale factor or the Hubble parameter has 
prevented any attempt to constrain the value of n against 
the observed data. Actually, only the model without 
matter has been investigated giving n G (—0.450, —0.370) 
or n G (1.366, 1.376) |30j | in clear disagreement with our 
estimate (|23H . However, such a comparison is meaning- 
less because of the presence of the matter term in the 
present case and the absence in the other one. 

Actually, using Ea. ((r3l) . it is possible to convert the 
estimate of n in a constraint on the present day value of 
the deceleration parameter. The best fit value for n thus 
translates into q$ = —0.33, while, combining Eas. i|15|) 
and $23$ , we get : 

( (-0.39,-0.27) at la 
q e\ . (26) 

[ (-0.43,-0.21) at 2a 
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FIG. 3: 1, 2, and 3<r confidence regions in the two dimensional 
parameter space (n,h). 



FIG. 4: 1, 2, and 3<r confidence regions in the two dimensional 
parameter space (n, £lb/&M)- 



While consistent with the picture of an accelerating uni- 
verse, our estimates for go disagree with other recent re- 
sults. Let us consider what is obtained for the ACDM 
model 3 . Using a flat geometry prior and fitting to the 
SNela Hubble diagram only, Riess et al. 0] have found 
n M = 0.29±£^ that gives 4 q = -0.56±0.07 that is 
not consistent with our estimate. Adding the data on 
the CMBR anisotropy and the power spectrum of SDSS 
galaxies, Tegmark et al. give Q M = 0.30±0.04 so that 
the estimated qo is in agreement with Riess et al. and 
hence in contrast with our value. A similar result has 
also been obtained by A04 only using the same f gas data 
we have considered here with a prior on h and fl^h 2 . For 
a flat ACDM model, their analysis gives fi M = 0.24±0.04 
and hence go = — 0.64±0.06 still in disagreement with our 
Eci. (|2t5[) . As a general remark, we notice that our mod- 
els turn out to be less accelerating (i.e., the predicted go 
is higher) than is observed for the standard concordance 
model. From a different point of view, lower values of 



3 We limit our attention to the ACDM model only since the cos- 
mological constant is the simplest and most efficient way to fit 
most of the astrophysical data 0. Moreover, the constraints on 
the equation of state parameter w = p/p a re s till consistent with 
the cosmological constant value w = —1 jjjjj. This conclusion 
is further strenghtened by the methods that aim at recovering 
the evolution of the dark energy density from the data in model 
independent way (see, e.g., |4(I and references therein). 

4 Hereafter, we will compute the error on qo propagating the max- 
imum la uncertainty on Qm- Although not statistically correct, 
this method gives a quick order of magnitude estimate of the 
error which is enough for our aims. 



qo correspond to higher n, i.e. to steeper power law La- 
grangians that are, however, disfavoured by the lower 
matter content of the corresponding best fit model. 

However, one could deem as unreliable the comparison 
among qo constraints obtained under different underlying 
cosmological models and look for model independent esti- 
mates of the deceleration parameter. For instance, Riess 
et al. have tried to constrain the deceleration parame- 
ter by using the simple ansatz q(z) = qo + {dq/dz) z= oz 
or resorting to a fourth order expansion of the scale fac- 
tor thus estimating also the jerk and snap parameters 
|4lj . While the (quite large) constraints on go shown in 
their Fig. 6 agree with our own in Eq. I|26(l , the vanish- 
ing of {dq/dz) z= o is clearly ruled out. It is interesting 
to notice, however, that a similar analysis performed in 
Ref. 01 expanding the scale factor up to the fifth order 
and using no priors at all gives different results. A glance 
at Fig. 2 in that paper shows that our ranges for go are 
indeed acceptable even if the best fit value quoted there 
(go = —0.76) is outside our 2cr interval. Moreover, Fig. 3 
of the same paper suggests that the jerk parameter is 
only weakly constrained and may be also consistent with 
a null value so that it is not possible to reject models 
with constant q{z). 

Actually, there is some evidence in favor of the model. 
First, the estimated Hubble constant is in good agree- 
ment with recent values quoted in literature. In the 
framework of the concordance model, a combined anal- 
ysis of the CMBR anisotropy spectrum measured by 
WMAP, the power spectrum of SDSS galaxies, the SNela 
Gold dataset, the dependence of the bias on luminosity 
and the Lya power spectrum lead Seljak et al. to finally 
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estimate h = 0.710±g;^f (at 99% CL) gg consistent 
with our range in Eq. l|24[) . Results in agreement with 
those of Seljak et al. (but with larger uncertainties) have 
also been obtained by applying the same method to less 
complete dataset and are not reported here for sake of 
shortness (see, e.g., 0,0 and references therein). It is 
even more appealing the agreement among our estimate 
of h and those coming from model independent methods. 
For instance, by combining different calibrated local dis- 
tance indicatores, the HST Key project finally furnish 
h = 0.72±0.08 [38| in quite good agreement with our re- 
sults. This conclusion is further strenghtened when com- 
paring to the results from time delays in lensed quasars 
|44j and Sunyaev- Zel'dovich effect in galaxy clusters [45| . 

Having constrained with the likelihood analysis both h 
and Qb/QMi we may derive J7m by using an independent 
estimate of fif,. Following Kirkman et al. [ 10{. we adopt : 

tt b h 2 = 0.0214±0.0020 
that, combined with Eas. (|24|l and Ea. (|25|l . gives: 

Q M = 0.28±0.02 (27) 

where the error has been roughly evaluated by propagat- 
ing the lcr uncertainties on h, flb/^M and fl h 2 approxi- 
mated as symmetric around the best fit values 5 . Eg. 127(1 
is in very good agreement with recent results. As already 
quoted above, using only the SNela Gold dataset, Riess 
et al. have found 0.29lg o3 for a flat ACDM model, while 
the analysis of Seljak et al. gives Om = 0.284°; ' 7 g 60 (at 
99% CL). Finally, fitting to the f gas data only with priors 
on both h and fish 2 , but not imposing the flatness con- 
dition ab initio, A04 estimates Qm — 0.245^q'q37, while 
including the CMB data, they get fi M = 0.26to.u4- All 



5 It is likely that this method underestimates the true error thus 
only giving an order of magnitude estimate. 



these results are in almost perfect agreement with our 
estimate of CIm which is indeed a remarkable success. 

Finally, we could use the estimated values of n, h and 
Qm and the age of the universe t to put constraints on 
the coupling constant (3 through Ea. l(14|l . However, this 
does not give us any useful information since we have no 
theoretical motivation that may suggest us what is the 
value of (3. On the other hand, the freedom we have in 
the choice of /3 leaves us open the possibility to find a R n 
model which fits both the SNela Hubble diagram and the 
fg as data and also predicts the right value of t$. 

B. f(R)=alogR 

Let us now discuss briefly the results for models with 
the logarithmic Lagrangian in Eq. I|18l) . With the follow- 
ing choice of the model parameters : 

VtMHla- 1 = 0.162 , h = 0.650 , Q b /Cl M = 0.184 (28) 

we get the best fit curves shown in Figs-ElandEl While 
both fits are indeed very good, it is interesting to note 
that the SNela Hubble diagram is now reproduced with 
great accuracy also for the two SNela with the highest 
redshift in contrast with what is observed for the power 
law Lagrangian models. This is likely a consequence of 
having this class of model a non constant deceleration 
parameter in agreement with what is suggested by Riess 
et al. (within the caveat noted above). 

Figs.0 and |H1 show the two dimensional projections of 
the 1, 2, and 3<r confidence regions on the subset pa- 
rameter space (SiM-ffoa -1 , h) and ((]Mi/ 2 a _1 ,J]i,/l]i}) 
respectively. It is clear that J]mHo q_1 anticorrelates 
with both h and fl^/fl m- From the projection on the 
(h, U /VIm) plane (not shown here), we see that these 
parameters are negatively correlated. Combining these 
plots, we may argue that the Hubble constant is posi- 
tively correlated with both Qm and a so that the an- 
ticorrelation with fij/i/oQ -1 is due to the degeneracy 
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FIG. 7: 1, 2, and 3cr confidence regions in the two dimensional FIG. 8: 1, 2, and 3<r confidence regions in the two dimensional 
parameter space (QmHqO. -1 , h). parameter space {QmHIoT 1 , fi(,/fiM). 



between h and a that turns out to be stronger than that 
between h and £Im- 

Let us now consider the constraints on the single pa- 
rameters. We get : 

( (0.148,0.174) at la 
tlMHla- 1 e { ; (29) 

[ (0.129,0.194) at 2a 



( (0.644,0.657) at la 
he I ; (30) 

[ (0.637,0.664) at 2a 



( (0.180,0.188) at la 

Q b /n AI e \ . (3i) 

[ (0.176,0.192) at 2a 

It is more useful to translate the constraint on the com- 
bined parameter VImHqOl~ 1 (whose physical meaning is 
not immediate) in a range for the present day value 
of the deceleration parameter. Using Eq.J2U, we get 
go = —0.55 as best fit value, while the confidence regions 
turn out to be : 

( (-0.56,-0.54) at la 
q e\ . (32) 

[ (-0.58,-0.52) at 2a 

Moreover, being q(z) no longer constant for this class 
of models, we may also estimate the transition redshift 



obtaining zt = 0.61 as best fit value and the following 
confidence regions : 

( (0.57,0.66) at la 
z T e{ . (33) 

[ (0.52,0.74) at 2a 

Even if the deceleration parameter is varying with the 
redshift z, our estimate of qo is still in disagreement with 
the estimates discussed in the previous subsection. As a 
general remark, we notice that, as for the class of models 
with power law Lagrangian, the estimated qo is higher 
than what is predicted by the best fit ACDM model. 
However, the disagreement is now less severe and, in- 
deed, a marginal agreement may be sometimes obtained 
by considering the 3a confidence regions. 

We may also compare the transition redshift that, 
for a flat ACDM model, is given by : zt — [2(1 — 
DmV^m] 1 / 3 — 1. Using, for instance, the estimate of 
given by Seljak et al., we get zt € (0.52,0.91) with 
zt = 0.71 as best fit in quite a good agreement with our 
Ea. (|33|) . Moreover, it is encouraging that our la confi- 
dence region has a non null overlap with that estimated 
by Riess et al., zt — 0.46±0.13, using the model inde- 
pendent parametrization of q(z) = qo + (dq/dz) z= QZ. 

Regarding the Hubble constant, the confidence regions 
for h are almost the same as those obtained for the power 
law Lagrangian case. Hence, we are still in agreement 
with previous results in literature. This is not very sur- 
prising since h is essentially determined by the fit to the 
low redshift SNela and, in this range, both Dl and Da 
are almost model independent. As a consequence, the 
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estimated h turns out to be the same whatever is the 
underlying cosmology and in agreement with what one 
should obtain by a linear fit to the z < 0.1 SNela data. 

From the constraints l|30|l and (|3 1|) and the value of 
£lbh 2 in Ref. we estimate: 

Q M = 0.27±0.03 (34) 

with the error evaluated as for that in Eq. I|27|) . This is in 
perfect agreement with both the result for the power law 
Lagrangian case and the other estimates quoted above. 

One could use Ea. (|34[l and the constraints on 
QmHqO.' 1 and h to narrow the range for the coupling 
parameter a. However, this does not give any useful in- 
formation since there is no way to theoretically predict 
the value of a. It is, on the contrary, more interesting 
to evaluate the present age of the universe using Eq. H20(l 
and the constraints in Eqs.J2Hl an d l|3t)(l . For the best 
fit model, it is to = 10.3 Gyr, while t ranges between 
10 and 11 Gyr for the parameters running in their la 
confidence regions. These values are too low when com- 
pared to the estimated to for the best fit ACDM model. 
For instance, the best fit vanilla model of Tegmark et al. 
predicts t = 13.54±g;|f Gyr more than 9a higher than 
our estimated upper value. Notice, however, that the dis- 
agreement is less severe (but still of high significance) if 
compared to to = 14.4~t\ 3 Gyr determined by Rebolo 
et al. [f| by fitting the ACDM model to the anisotropy 
spectrum measured by WMAP and VSA and to the clus- 
tering properties of 2dFGRS galaxies. However, even if 
in agreement with those obtained by completely different 
methods, these estimates are model dependent. Actually, 
our predicted to is not unreasonably low if we consider 
that globular clusters data lead to to = 12.61^6 Gyr 
|47| . while a lower limit to > 12.5±3.5 Gyr is obtained 
by nucleochronology . Considering the 2a confidence 
regions for the parameters f^/ii^a -1 and h, it is there- 
fore possible to find models that are able to successfully 
fit the astrophysical data we are considering (even if they 
are not the preferred ones) and also predict a present age 
of the universe that is not in disagreement with cosmol- 
ogy independent estimates of to- 

V. DISCUSSION AND CONCLUSIONS 

Assuming that the Palatini (first order) approach is 
the correct way to treat /(i?) theories, we have investi- 
gated the viability of two different class of cosmological 
models corresponding to two popular choices of f(R), 
namely a power law in the scalar curvature and a loga- 
rithmic function of R. The expansion rate H — a/a may 
be analytically expressed as a function of the rcdshift z 
for both classes of models so that it is possible to contrast 
the model predictions against the observations. In par- 
ticular, we have used the SNela Hubble diagram and the 
data on the gas mass fraction in relaxed galaxy clusters 
to investigate the viability of each class as dark energy 



alternative and to constrain their parameters. The main 
results are sketched below. 

1. Both classes of models provide very good fits to 
the data even if the choice f(R) = alni? leads 
to a Hubble diagram that is in better agreement 
with the trend shown by the highest redshift SNela. 
However, the paucity of the data does not allow us 
to eventually prefer one model to the other. 

2. Eas. (|22|l and (|28fl give the best fit parameters for 
the power law and logarithmic Lagrangian mod- 
els respectively. The confidence regions have been 
determined from the marginalized likelihoods and 
are reported in Eas. 12311 - 125|) for the models with 
f(R) = l3R n and in Eqs.® - (EB for those with 
f(R) = alni?. To better compare the model pre- 
dictions with previous results in literature, we have 
evaluated the present day deceleration parameter 
go and the matter density parameter Qm (assuming 
the estimate of Qi,h 2 in Ref. 46]). For both classes 
of models, qo turns out to be higher than what is 
predicted for the concordance ACDM model, i.e. 
f(R) theories lead to models that accelerate less 
than what is observed. This result is however some- 
what weakened by comparing with model indepen- 
dent estimates of qo even if these latter may be 
affected by systematic errors. As far as the mat- 
ter content is concerned, for both classes of models 
f2jvf is in very good agreement with what is inferred 
from galaxy clusters and estimated by fitting the 
ACDM model to the available astrophysical data. 

3. To ameliorate the agreement with the observed qo, 
one should increase the value of n for the mod- 
els with power law Lagrangians or decrease that of 
^mHqC -1 for models with /(i?) = alni?. In this 
case, a good fit to the data may still be obtained 
provided that both h and Qb/^M are increased. 
While higher values of h could still be compatible 
with the local estimates of the Hubble constant, in- 
creasing ilb/flM leads to lower values of FIm ■ Actu- 
ally, the very good agreement among the estimated 
Qm in Eas. (|27|l and (|34|l and the results in litera- 
ture is a strong evidence against this choice. There- 
fore, we conclude that it is not possible to recover 
the same value of qo in the concordance model by 
using power law or logarithmic Lagrangians. 

4. A model independent estimate of the present day 
age of the universe to allows one to break the mat- 
ter/geometry degeneracy inherent in f(R) theories 
recovering the value of the coupling constant. For 
power law Lagrangians, this is indeed the only way 
to determine (3 thus offering the possibility to al- 
ways recover a model that both fits the SNela Hub- 
ble diagram and the data on the gas mass fraction 
in relaxed galaxy clusters and also has the correct 
age. On the other hand, to is an independent check 
for models with logarithmic Lagrangian since, in 
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this case, it may be evaluated as a function of the 
two parameters VLmH^oi^ 1 and h and compared 
with previous results in literature. It turns out that 
the predicted to is lower than the value estimated 
for the ACDM model and only marginally consis- 
tent with what is inferred from globular clusters 
and nucleochronology. 

The results summarized above may pave the way to 
the solution of an intriguing dilemma : is Einsteinian gen- 
eral relativity the correct theory of gravity? If yes, then 
dark energy is absolutely needed to explain the accel- 
erated expansion of the universe and hence all the the- 
oretical efforts of cosmologists have to be dedicated to 
understanding its nature. On the contrary, if f(R) theo- 
ries are indeed able to explain the accelerated expansion, 
then it is time to investigate in more detail what is the 
right choice for the function f(R) and how the variation 
has to be performed (higher order metric or first order 
Palatini approach). 

From the observational point of view we have adopted 
here, there are no strong evidences against models with 
power law or logarithmic Lagrangians in the framework 
of the Palatini approach. On the contrary, we have 
seen that both classes of models successfully fit the data 
with values of the Hubble constant and matter content 
in good agreement with some model independent esti- 
mates. However, there are some hints that could lead 
to reject both choices for f(R). Models with power law 
Lagrangians have a constant q(z) so that they are always 
accelerating. This is not consistent with the (tentatively) 
observed transition from acceleration to deceleration at 
zt — 0.5. Moreover, a constant q{z) could give rise to 
problems with nucleosynthesis and structure formation. 
On the contrary, models with a logarithmic Lagrangian 
are not affected by such problems and indeed they pre- 
dicts a transition redshift which is in good agreement 
with the estimates for the ACDM model. On the other 
hand, these models turn out to be too young, i.e. to is 
lower than what is expected. 

Actually, a more general remark is in oder here. Let us 
suppose we have found that a given choice for f(R) leads 
to models that are in agreement with the data so that 
we should conclude that this class of models correctly 
describe the present day universe. What about the early 
universe? One could expect that the functional expres- 
sion of f(R) is not changing during the evolution of the 
universe, even if R may evolve with cosmic time. If this 
were the case, then the correct choice for f(R) should 
be the one that leads to models that are not only able to 
reproduce the phenomenology we observe today, but also 
give rise to an inflationary period in the early universe. 
Therefore, we should reject logarithmic Lagrangians since 
it is well known they do not predict any inflationary pe- 
riod. On the other hand, the choice f(R) = (3R n is able 
to explain inflation provided one sets n = 2, not too far 
from our estimate in Ea. i|23|) . From this point of view, 
it is worth noticing that the astrophysical data we have 
considered probe only the present day universe, while to 



depends on the full evolutionary history. Indeed, loga- 
rithmic Lagrangians fail to reproduce the correct to in 
the same way as they fail to give rise to inflation, while 
both inflation and to are correctly predicted by models 
with power law f(R). This may argue in favour of this 
choice for f(R), but actually there is no reason to ex- 
clude the possibility that also the functional expression 
of f(R) changes with time so that neither class of models 
may be definitively rejected or deemed as the correct one 
from this point of view. 

Summarizing, the likelihood analysis presented here al- 
lows us to conclude that the Palatini approach to f(R) 
theories leads to models that are able to reproduce both 
the SNela Hubble diagram and the data on the gas mass 
fraction in galaxy clusters. From an observational point 
of view, this means that both power law and logarithmic 
f(R) are viable candidates to explain the observed accel- 
erated expansion without the need of any kind of dark 
energy. However, open questions are still on the ground. 

First, we have not yet been able to discriminate be- 
tween the two classes of models. Theoretical considera- 
tions and some hints from the age of the universe could 
argue in favour of the power law f(R), while the observed 
transition from acceleration to deceleration in the past 
disfavors this choice. To solve this issue, one has to re- 
sort to high redshift probes such as the CMBR anisotropy 
spectrum. While the data are of superb quality, the un- 
derlying theory is still to be developed so that fitting the 
CMBR anisotropy temperature and polarization spectra 
with f(R) theories will be quite a demanding task. 

Second, we have only considered two physically moti- 
vated and popular choices for f(R). Several other mod- 
els are possible and are worth of being tested against the 
data. In particular the R\nR Lagrangian which is re- 
lated to the Straobinsky inflationary model [49( and to 
the limit R n -> R for n ~ 1 being 

= RR e — R ( e e ln R ) ~ R + e R In R + 0(e 2 ) . (35) 

However, rather than being confused by a plethora of 
successfull models, it is desiderable to develop a method 
that allows to directly reconstruct f(R) from the data 
with as less as possible aprioristic assumptions. This will 
be the subject of a forthcoming paper |51| . 

Last but not least, whether the Palatini approach is 
indeed the correct method to treat f(R) theories or the 
metric approach should be preferred is still an unsolved 
problem. We have shown here that the Palatini approach 
is not rejected by the data, but a similar analysis for the 
same models considered in the framework of the metric 
approach is still lacking. However, it is worth noticing 
that even this test will not be conclusive. Let us con- 
sider, for instance, two choices fi(R) and fi{R) and let 
us suppose that fi(R) fit the data if considered in the 
framework of the metric approach, but not if the Palatini 
approach is used. Let us further assume that the opposite 
holds for f2(R). From an observational point of view, it 
is impossible to select between fi(R) and f2(R)- Hence, 
observations could never suggest what is the correct way 
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of performing the variation of a f(R) Lagrangian. The 
answer to this question is outside the possibilities of an 
astronomer and lies fully in the field of a theoretician. 

As a final comment, we would like to stress the need 
for synergy between theory and observations. While it is 
possible to build a physically motivated and mathemati- 
cally elegant theory, it is not so easy to fit the significant 
amount of astrophysical data now available. Since the 
words observational and cosmology may today be joined 
together in a single meaningful term (observational cos- 
mology), it is time to look at every theoretician's proposal 
from an observational point of view before drawing any 
conclusion about the validity of a whatever model. Even 



if not always conclusive, in our opinion, this is the only 
way to shed light on the dark side of the universe. 
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